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The existence of Dark Matter (DM) in the form of Strongly Interacting Massive 
Particles (SIMPs) may be motivated by astrophysical observations that challenge 
the classical Cold DM scenario. Other observations greatly constrain, but do not 
completely exclude, the SIMP alternative. The signature of SIMPs at the LHC may 
consist of neutral, hadron-like, trackless jets produced in pairs. We show that the 
absence of charged content can provide a very efficient tool to suppress dijet back¬ 
grounds at the LHC, thus enhancing the sensitivity to a potential SIMP signal. We 
illustrate this using a simplified SIMP model and present a detailed feasibility study 
based on simulations, including a dedicated detector response parametrization. We 
evaluate the expected sensitivity to various signal scenarios and tentatively consider 
the exclusion limits on the SIMP elastic cross section with nucleons. 
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I. INTRODUCTION 

There is ample evidence of the existence of Dark Matter (DM) in the Universe. First 
invoked to explain a puzzle in clusters of galaxies [T], the case for DM, based on obser¬ 
vations from galactic up to the largest scales, is now very strong (see for instance BE]). 
In particular, precision measurements of the cosmic microwave background anisotropies 
imply that about 80% of the matter content of the Universe is made of DM [H E], De¬ 
spite this, the precise nature of DM remains a mystery. The most-studied hypothesis 
assumes that DM is made of weakly interacting massive particles or WIMPs. This rests 
on the observation that a stable massive particle with annihilation cross section of or¬ 
der 1 pb -characteristic of weak interactions- could have a relic abundance that agrees 
with cosmological measurements. A popular WIMP is the neutralino, the typical DM 
candidate of supersymmetric extensions of the Standard Model (SM) of particle physics 
[B], but many alternative WIMP candidates have been proposed [7], which are actively 
being searched for by experiments, including those at the LHC. Yet, as appealing as the 
WIMP scenario may be, it remains important to study other possibilities. Amazingly 
enough, a strongly interacting massive particle (SIMP), i.e. a particle with strong inter¬ 
actions with ordinary baryons, is not yet fully excluded [S] (see also [B] and specially BDl 
for a more recent appraisal of existing constraints). While a SIMP scenario may seem 
exotic at hrst sight, it may be motivated by the long lasting interest for DM particles 
with strong self-interactions, going back to the seminal work of Spergel and Steinhardt 
m Indeed, self-interacting dark matter (SIDM) particle^ with very large cross section, 
^xx/'^x ~ 10“^^ cm^/GeV, may help addressing astrophysical observations that present 
a challenge for the cold dark matter (CDM) paradigm, like the missing satellites or core¬ 
cusp problems [IMS]. While it is possible to build scenarios with a strongly interacting 
hidden sector weakly coupled to the SM particles (see e.g. dl). it is perhaps as natural 
to consider SIDM particles with strong interactions with ordinary matter. The latter 
possibility is much more constrained but, again, it is not fully excluded [Sj. This will be 
our main motivation to take the SIMP hypothesis seriously. In any event, the absence of 
a clear signal in recent WIMP searches may be a further motivation to look away from 
the lamppost and to explore a bit more this exotic possibility. Here we focus in particular 
on the possible search of a SIMP candidate at the LHC. 

In particular we consider the possibility of observing events with no (or little) signal 
in the tracking systems and electromagnetic calorimeters and only (or essentially) energy 
deposition in the hadron calorimeters of the detectors, akin to the signature of neutrons 
or K^. As far as we know, this possibility has been hrst put forward in [18], albeit in 
the framework of the Tevatron, and more recently in [12]. We will follow the simplihed 
model approach of the latter work but our feasibility study goes quite a few steps further. 
In particular, we develop the charged contents as a discriminator to suppress dijet back¬ 
grounds at LHC, thus enhancing the sensitivity to a potential SIMP signal. We present a 
feasibility study based on simulations, including a Delphes 123 description of a typical 
LHC detector, and evaluate the expected sensitivity to various signal scenarios. We h- 
nally tentatively map the expected sensitivity onto the elastic cross section for scattering 
of SIMPs on nucleons and compare the resulting exclusion limit (i.e. assuming no signal is 
seen at the LHC) to other constraints, mostly astrophysical but also from direct detection 

^ Unfortunately, strongly interacting and self-interacting share the same acronym, SI, so that in the 
literature SIDM may stand for strongly interacting dark matter or SIMP for self-interacting massive 
particles (see e.g. m)- 
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DM searches. 

The plan of this article is as follows. Firstly, we briefly motivate the models we con¬ 
sider in the light of the existing constraints (Section]^. Next we describe the expected 
signature at the LHC, and the dedicated simulation of the signal and background sam¬ 
ples (Sections III A and IIIB), and detail the strategy of our analysis (Section IIIC), 


before giving the resulting expected sensitivity on various signal scenarios (Section HID). 
In the last section, we tentatively re-express this in terms of an exclusion limit on the 
SIMP-nucleon cross section (Section |lv]) and hnally draw our conclusions (Section [V|). 

Before to go on, we mention that other scenarios with a strongly coupled hidden sector 
that may lead to LHC signals somewhat analogous to trackless jets have been proposed 
recently, see [2I1E2]. If necessary, those may provide a further motivation for the kind of 
experimental study we discuss in the following sections. 


II. SIMPLIFIED SIMP MODELS 

In this section we motivate the SIMP models that we consider. Taking into account 
the existing constraints, the possibilities boil down to a few options m, at least provided 
we focus on simple models. By this we mean assuming that the SIMP, which may be 
composite, can nevertheless be treated as an elementary particle on all energy scales we 
deal with. We will furthermore suppose that the SIMP particles interact with SM particles 
(here quarks) through a mediator elementary particle. In short, we apply the philosophy 
of simplihed models used for WIMP searches at colliders ^31 to SIMP phenomenology 
m- Given the exploratory character of the search for DM it may be reasonable to begin 
with such simple assumptions. Whether the conclusions that can be drawn using such 
a framework are generic is another question. The dynamics of strong interactions of the 
ordinary hadrons is very complex. The same is probably true of realistic SIMP scenarios, 
see e.g. 12a. Yet, the phenomenological approach we advocate allows to go signihcantly 
beyond the approximation usually adopted to describe SIMP interactions. In particular, it 
allows to confront SIMP properties advocated to solve astrophysical issues, with possible 
signatures at high energies. In the rest of this section, we summarize the properties of 
our simplihed SIMP models, following essentially the arguments of [T^ . 

One of the strongest constraints on a SIMP as DM is set by searches for heavy isotopes, 
in particular heavy water, which put limits on the formation of bound states between 
SIMPs and nucleons. Assuming that the SIMP is the dominant form of DM, a particle 
lighter than ~ 10 TeV that can form a bound is excluded [S] (see also Ea). This constraint 
is avoided if the SIMP has purely repulsive SIMP-nucleon interactions, which may be 
achieved if it interacts with SM particles through a scalar or vector mediator with opposite 
sign couplings. In the case of a vector mediator we should be concerned with the fact 
that vector mediators couple to DM antiparticles with an opposite charge [12]. This 
is avoided if there are no DM antiparticles around, that is if the abundance of DM is 
asymmetric. This could actually be the case. First, a symmetric SIMP candidate can 
only be a sub-dominant component of DM if its abundance is set by thermal freeze-out. 
Conversely, if most of DM is made of SIMPs, then its abundance is determined either 
by an asymmetry or through a non-thermal mechanism. Second, measurements of the 
Earth heat flow set strong constraints on SIMP properties. For cross sections that are 
characteristic of a SIMP, DM could be efficiently captured and accreted in the core of 
the Earth. Annihilating SIMPs would provide a substantial source of heat, a constraint 
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that does not apply to asymmetric candidates [in]|^ In the sequel we will take seriously 
the possibility that DM is made dominantly of a SIMP, if anything because it may be of 
interest to compare the LHC reach with the expectations from other DM searches, like 
direct detection ones for example, which implicitly rely on this assumption. 

The interaction Lagrangian of the models we will consider is then simply (see also m) 


-Qx^t^XX - 


qq 


(scalar mediator) 


An* = <! (1) 

~9x4>ii Xl^X ~ 9q<PtJ. qi^q (vector mediator) 

with gx9q,9x9q < 0 to avoid the formation of bound states. Further parameters of the 
model are the mass of the messenger, and of the SIMP, m^. The model has thus 
at least 4 free parameters. For LHC phenomenology, only the product of the two cou¬ 
plings appears, but astrophysics constrains both DM self-interaction and interactions with 
ordinary matter. Furthermore, although we have in mind the fact that SIMPs shoulcQ 
have flavour dependent couplings, for simplicity we assume in this study that they have 
a universal coupling to quarks. 

Introducing new strong interactions between quarks, and thus nucleons, is not harmless. 
In the sequel we will consider a rather light mediator, ~ 1 GeV, as in [12] . While the 
precise value of this mass may seem a priori of little importance for SIMP production at 
the LHC, we will argue that it is crucial to assess their detection in the hadron calorimeters 
as well as for the comparison with other DM searches. In order to keep small the impact 
of the new interaction on the nuclear potential, we will furthermore assume, again as in 
ESI, that the mediators do not modify nuclear potentials by more than O{10%), so that 
9xN ^ ^■^9nN 3 for m</, ~ ttIt,, where g^N ~ 13 is the effective pseudoscalar pion-nucleon 

coupling [22|, or, for 1 GeV, g^N < O.SgpN ~ 6, where gpN 18 |3U] is the vector p 

meson coupling to the nucleon. This will be relevant in Section IV when we will discuss 
the comparison between high energy and low energy constraints. 

There are further constraints on the interaction between DM and ordinary matter, and 
between DM particles themselves that we should take into account. First, one may look 
for elastic scattering of SIMPs with nuclei in direct detection experiments. However, as 
a SIMP interacts strongly in the Earth (or even in the atmosphere), it cannot reach the 
underground direct detection detectors [21] Instead, such high interaction cross sections 
are constrained by space or airborne experiments like RSS, a balloon-based detector [22] 
and XQC, a sounding rocket X-ray experiment [21]. There are also constraints from 
primordial nucleosynthesis and cosmic rays [23 ES]- these constraints have been 
extensively reviewed in EDI, to which we refer for more details since we have nothing 
specihcally new to say about them. Finally, there are strong constraints on interactions 
between DM and baryons from observations of Cosmic Microwave Background Radiation 
(CMBR) anisotropies and large structure, including from Lyman-a data [23 EB]- In 
particular the constraints reported in [2E] are relatively new and somewhat stronger than 
previously thought. We will discuss them further in section IV where we also discuss the 
cosmological constraints on DM (self)interactions. 


^ Notice however that, in an asymmetric scenario, light scalar DM particles can lead to black hole 
formation if they are trapped inside neutron stars BZIISE]. This is valid for light asymmetric bosonic 
DM candidates. For simplicity, we only consider fermionic DM candidates. 

^ Light SIMPs significantly coupled to b or c quarks is probably constrained by B and D meson phe- 
nonomenology. 

^ Some interesting implications for direct detection experiments of significant energy loss of DM in the 
Earth are discussed in [5^ . 
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The possible signatures of SIMPs at colliders have been much less studied [TSl ITU] , 
The production at LHC of WIMPs is currently studied through missing transverse energy 
signals, typically jets or photons with large missing momentum Such constraints 

also apply to SIMPs, provided their interaction with baryons is less than that characteris¬ 
tic of hadrons, so that the DM particles produced do not deposit substantial energy in the 
calorimeters of the detectors. In the present work, we discuss further the complementary 
possibility of observing trackless jets from DM interactions. 


III. SIMP OBSERVABILITY AT THE LHC 

The potentially high interaction cross section of the SIMPs with nuclear matter leads 
to a particular signature at the LHC. Indeed, the neutral SIMP leaves no track in a 
tracking detector, but will give rise to a highly-energetic hadronic shower. Since a pair 
of SIMPs will be dominantly produced back-to-back in the transverse plane, this leads to 
the peculiar observable signature of a dijet pair without any tracks. Here, we study the 
feasibility of detecting this signature at the LHC, extending on the previous work in [TU] . 


A. Event generation 

The interaction Lagrangian Q was implemented in FeynRules 2.0 [H] for the scalar 
mediator case and subsequently interfaced to Madgraph 5 |1S] to generate the pair 
production of SIMPs in proton-proton collisions. The center-of-mass energy was chosen 
to be ^/s = 8 TeV, corresponding to the energy at which the LHC delivered collisions in 
2012. The benchmark for our simulation is defined by the couplings = —1, Qq = 1 and 
the mediator mass to = 1 GeV (the precise value of plays no role at this level of 
the discussion). We consider a stable dark-matter particle y, and generated events for 
various masses = 1, 10, 100, 200, 400, 700, and 1000 GeV. 

In Table m the SIMP production cross section is shown for each considered SIMP mass 
rriy^, along with the number of events expected for an integrated luminosity of 20fb~^, 
which more or less corresponds to the dataset recorded by the LHC experiments in 2012. 
These values were obtained using |p(x)| < 2.5 and Pt{x) > 250 GeV pres-election re¬ 
quirements. Along with the case of a scalar mediator, the consistently larger values for 
a vector mediator are reported as well. For this LHC feasibility study, we have focused 
our simulations on the scalar case, but the results can be directly translated to a vector 
mediator by straightforward scaling of the reported cross sections. 

At these high cross sections of dijet-like signal events, we consider QCD jet production 
as the main background process to take into account. Also this background was generated 
with the MadGraph program, where we produced two samples in the ]500,1000] and 
jlOOO, inf[GeV bins of = Z]partons |PT,i|, the scalar sum of the transverse momenta of 
the outgoing partons. The corresponding cross sections are respectively 8426 and 204 pb. 
With the above settings, the possible additional production of dijet events through the 
mediator 0 was verihed to be at the percent level with respect to the QGD dijet production, 
and is further ignored. 

Events from both signal and background samples are subsequently processed with 
PythiaS [in], using the CTEQ6L1 El parton distribution functions, in order to embed 
the hard interactions in full proton collisions, including the description of the parton 
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Scalar 

Vector 

[GeV] 

[Pb] 

Events/20fb ^ 

f^xx [Pb] 

Events/20fb ^ 

1 

2.80 

56040 

3.18 

63645 

10 

2.79 

55800 

3.17 

63400 

100 

1.88 

37620 

2.49 

49826 

200 

0.728 

14558 

1.31 

26196 

400 

0.0769 

1539 

0.229 

4583 

700 

0.00363 

73 

0.0167 

336 

1000 

0.000239 

5 

0.00147 

31 


TABLE 1. Production cross section for each SIMP mass, and number of events corresponding 
to an integrated luminosity of 20fb“^, after \r]{x)\ < 2.5 and Pt{x) > 250GeV pres-election 
requirements. 

shower, hadronization, and underlying event. 


B. Detector simulation 

In order to study the observability of a SIMP signal at the LHC, we use the parametrized 
detector simulation package Delphes [20] • Using this program, we simulate the response 
of a typical LHC detector to the generated signal and background collisions. We use 
Delphes in the standard CMS |1H] conhguration, which implements spatial and energy 
resolution functions for each sub-detector — inner tracker, electromagnetic and hadronic 
calorimeter, and muon spectrometer — and adds high-level event reconstruction, like jet 
clustering using Fast Jet [1H]. Some parts of this event reconstruction have particular 
relevance for this analysis. 

First, the track reconstruction plays an important role in establishing jets to arise from 
only neutral particles. In Delphes, tracks are built from generated particles, to which a 
realistic inefficiency function is applied. Hence, only genuine tracks are being simulated. 
This is adequate for our purpose, since we expect our QCD dijet background to arise from 
a combination of tracking inefficiency and of fluctuations in the jet fragmentation, leading 
to a genuinely small charged jet content. 

Further, an adequate simulation of the jet response is needed. This is obtained in 
Delphes through the smearing of jet energy measurements with realistic jet resolution 
functions, rather than using a detailed simulation of the calorimeter response. For the 
SIMPs, though, we implemented a more elaborate approach. The interaction of the SIMP 
X with the detector’s calorimeters can be described as an interaction of a hadron with the 
nuclei of the detector material, but comes with a potentially different nucleon interaction 
cross section We model this interaction in a way that allows to easily simulate a 

change of ct^at, for instance leading to incomplete containment of the deposited energy in 
the calorimeters. 

The position of the first yiV collision serves as the starting point of an ensuing hadronic 

ln(2) ^ 

shower. This initial position is distributed exponentially as e ^ where x is the 

calorimeter depth in interaction lengths A/ and ctqcd is the standard hadronic interaction 
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cross section of ~ 40 mb, which is taken constant in the energy range of the collisions 
of interest. We study the dependence of the search sensitivity on by considering 
c’‘xA^/‘^QCD = 1, where full shower containment is expected, and ct^at/ctqcd = 0.1, where 
late shower development will lead to calorimeter energy leakage. As will be detailed later, 
this lower value allows to gauge the dependence of the analysis on ct^n- A larger value of 
was not explicitly studied, but can up to some level still be assumed to lead to full 
registration of the energy of the induced early showers. 

After the hrst collision, the shower will develop as a mixture of subsequent collisions 
of the SIMP with standard hadronic collisions of particles created in the shower. For 
simplicity, we assume standard longitudinal hadronic shower development. This is ex¬ 
pected to be a good approximation, since the energy loss of the SIMP is largest in the 
beginning of the shower. We model the longitudinal energy prohle of the shower according 
to |Sni, with parameters estimated for iron absorber 1^ . The total deposit of energy 
in the calorimeters is then calculated as an integration of the longitudinal energy prohle, 
from the starting position of the shower up to the rear face of the calorimeter. For this, 
we assume a calorimeter depth of 1.7A/ and 9 A/ for our electromagnetic and hadronic 
calorimeters respectively, uniform as a function of polar angle of the incident particle. 

Finally, due to the high instantaneous luminosities reached by the LHC, additional 
proton-proton interactions, called pile-up, may occur simultaneously in each LHC bunch 
crossing. These pile-up collisions are also simulated by Delphes, where we set the average 
number of pile-up interactions to 21, Poisson distributed, and with a resolution on the 
z coordinate of the interaction vertices cTz = 0.1mm. This conhguration approximates 
sufficiently the reality of the data-taking in the 2012 LHC run. 

For events containing pile-up interactions, we apply the Delphes pile-up subtraction 
algorithm in order to remove the additional particles and their energy deposits from the 
reconstruction. This is, on the one hand, important to avoid charged particles from 
the pileup interactions that overlap with SIMP jets to induce inefficiencies in the signal 
reconstruction. It was checked that this does not happen at an appreciable level. On 
the other hand, the energy deposits from pile-up collisions may bias jet energy estimates, 
affecting the analysis selection. Also this effect was checked, and the pile-up subtraction 
was found to remove any such signihcant pile-up dependence. 


C. Analysis 

The SIMP signal considered in this paper would manifest itself in the detector as a pair 
of trackless jets, back-to-back in the plane transverse to the beamline. The background 
to such a signal is expected to be dominated by standard QCD dijet production, where 
an interplay of rare jet fragmentation and tracking inefficiency leads to a very low number 
of tracks in the jet, and a very low fraction of jet energy emitted in charged particles. To 
quantify the latter, we use as an observable the so-called charged energy fraction (CHEF), 
dehned as the ratio I]iPT,i/PTjet, where the sum runs over all tracks with transverse 
momenta pT,i, associated to the jet with transverse momentum PTjet- 

The selection of the signal events with high-px back-to-back jets proceeds as follows. 
At least two jets are required with px > 350 GeV and |pjet| < 2.0, with an azimuthal 
separation | A0jj | > 2.0 on the two leading jets. The px requirement on the jets is driven by 
a typical dijet trigger requirement for 2012 LHC conditions, while the r] restriction ensures 
the tracks of the jets to be well contained within the tracking detector. In Figure [T| the px 
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FIG. 1. The jet pT of the two leading jets after the selection cuts described in the text, excluding 
the jet pt cut, for cT^iv/(TQCD = 1 (left) and /cjqcd = 0.1 (right). In all cases, a generator- 
level pt pres-election at 250GeV is already applied. 


of both leading jets is shown after the above cuts, excluding the pt cut, for all considered 
signals and for the background, for both cases o'xn/o'qc'd = 1 and o'xn/o'qcd = 0.1. The 
similarity of the shape of signal and background distributions at high px is a manifestation 
of the lack of possible discrimination of signal and background on a purely kinematic basis. 

The number of tracks or the CHEF of a jet provide strong handles to suppress the 
background to the SIMP signal. In the remainder, the CHEF is used, as it was found 
to have a better performance than the number of tracks. The distribution of the CHEF 
observable is shown in Figure]^ (left) for the leading jet, comparing signal to background, 
after the kinematic cuts above. In Figure (right), the background suppression is com¬ 
pared to the signal selection efficiency, for various cuts on the CHEF of both the two 
highest-pT jets. A background rejection with a factor of 10000 can be achieved with a 
signal selection efficiency above 95%. 

While the CHEF is the main handle to separate a SIMP signature from QCD jets, 
other differences between signal and background exist. These arise due to the nature 
of SIMP showers being different than the QCD background jets, the latter arising from 
quark or gluon fragmentation. The SIMP, being a single particle, is expected to induce a 
narrower hadronic shower, potentially with a different longitudinal shower development, 
leading to differences in the shower shape, the jets’ electromagnetic fraction, or giving rise 
to increased punch through of particles from the hadronic showers in the muon system. 
Although exploiting such effects can significantly enhance further the analysis sensitiv¬ 
ity |S3], it requires more detailed simulations, beyond the scope of this study. 

The cut on the CHEF of the two leading jets is chosen to ensure a large signal-to- 
background ratio over a wide range of SIMP masses. Contrary to many searches at the 
LHC, not much analysis optimisation is possible, nor needed. On the kinematic side, 
the cuts are fully defined by trigger and detector acceptance, as described above. What 
concerns the CHEF cut, it is mostly defined by the available integrated luminosity, since 
the production cross section of the SIMPs cannot vary much without losing the signal, 
despite a large signal-to-background cross section ratio. Indeed, as the SIMP inelastic 
interaction and related production cross sections decrease, the signal will quickly transition 
from a visible pair of trackless jets into missing energy, and the analysis presented in this 
paper unavoidably loses its sensitivity. In order to provide a way to gauge this effect. 
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FIG. 2. Left: charged energy fraction of the highest-px jet; after kinematic selection requirements 
as described in the text. Right: Background suppression factor versus signal selection efficiency 
for various cuts on the CHEF, labeled on the graph. 


we assessed the search sensitivity on the samples simulated with a factor 10 reduction in 
the inelastic interaction cross section which is proportional to Because of the 

same proportionality, this leads to an equal reduction by a factor of ten of the production 
cross sections mentioned in Table This latter reduction is further always implied when 
considering a^N/(T qcd = 0.1. 


D. Results 


In Table we snmmarize the nnmber of events after the various analysis cuts, for all 
signal and background samples, assuming ct^n/o'qc'd = 1- Table III contains the results 
for cr^Ar/cjQCD = 0.1. To overcome the statistical limitation from the hnite simulated back¬ 
ground samples, we calculated the number of QCD multijet background events passing all 
cuts by assuming uncorrelated efficiencies of the CHEF requirements on the two leading 
jets. We checked this assumption to be conservative, in the sense that for loose cuts on the 
jets’ CHEF, the assumption of uncorrelated cuts leads to a similar or smaller background 
reduction compared to a direct estimation from applying the CHEF requirement on both 
jets in the simulated events. 

The results show readily that the considered SIMP signal can be discovered at the 
LHC for low SIMP mass. At high mass, the sensitivity fades out as the signal production 
cross section drops. We calculate the signal signihcance as a function of SIMP mass using 
the Z-value given by [M] 

Z = V2eTr\l - 2p) (2) 


with p-value 


p = B 


1 + r 


,s + b,Tb + l 


(3) 


where B is the incomplete beta function dehned in [5H], s is the number of signal events, 
b is the number of background events, and r is the number of background events divided 
by the background uncertainty squared. For this study, we do not provide estimates of 
the expected precision on the background prediction, but rather consider three levels of 
uncertainty, 20%, 50%, and 100%, as such providing a range in which the impact of the 
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Sample 

Total 

Selection 

Gharged Energy Eraction 

0.1 

0.05 

0.04 

0.03 

QCD HT- 5 OOT 0 IOOO 

168520000 

4992040 

378 

45 

21 

11 

QCD HT-lOOOToInf 

4080000 

2691870 

144 

16 

8 

3 

SIMP (1 GeV) 

56040 

9076 

8693 

8553 

8504 

8436 

SIMP (10 GeV) 

55800 

9060 

8674 

8531 

8479 

8413 

SIMP (100 GeV) 

37620 

7252 

6950 

6841 

6804 

6748 

SIMP (200 GeV) 

14558 

4072 

3928 

3871 

3848 

3821 

SIMP (400 GeV) 

1539 

771 

749 

739 

735 

730 

SIMP (700 GeV) 

73 

52 

51 

50 

50 

49 

SIMP (1000 GeV) 

5 

4 

4 

4 

4 

4 


TABLE 11. The number of events after the selection cuts and several cuts on the charged energy 
fraction, for the 2 background samples and the 7 different signal samples with cr^Ar/cjQCD = 1- 


Sample 

Total 

Selection 

Gharged Energy Eraction 

0.1 

0.05 

0.04 

0.03 

QGD HT- 5 OOT 0 IOOO 

168520000 

4992040 

378 

45 

21 

11 

QGD HT-lOOOToInf 

4080000 

2691870 

144 

16 

8 

3 

SIMP (1 GeV) 

5604 

146 

134 

131 

131 

130 

SIMP (10 GeV) 

5580 

144 

132 

130 

129 

128 

SIMP (100 GeV) 

3762 

117 

107 

105 

105 

104 

SIMP (200 GeV) 

1456 

68 

63 

62 

62 

61 

SIMP (400 GeV) 

154 

14 

13 

13 

13 

13 

SIMP (700 GeV) 

7 

1 

1 

1 

1 

1 

SIMP (1000 GeV) 

1 

0 

0 

0 

0 

0 


TABLE III. The number of events after the selection cuts and several cuts on the charged energy 
fraction, for the 2 background samples and the 7 different signal samples with (T^^tv/ctqcd = 0.1. 


background uncertainties is demonstrated. In Tables |IV| and [V] the Z-value is given for the 
considered background uncertainties and several CHEF cuts, respectively for the SIMP 
samples with cr^Ar/c’‘QCD = 1 and cr^Ar/cTQCD = 0.1. 

We can conclude from Table [TV] that a discovery can be made through the observation 
of a 5cr excess for SIMP masses up to 400 GeV, for a CHEF cut of 4% or tighter, for all 
considered background uncertainties, and for QCD-level SIMP interaction cross section 
This result for a scalar mediator holds as well in the case of a vector mediator, for 
which a small cross section increase needs to be accounted for, as reported in Table [1} 
For the o-^^joQCD = 0.1 case, with a background uncertainty of 20% or 50%, a CHEF 
cut of 3% is needed in order to reach discovery up to = 100 GeV. If the background 
uncertainty amounts to 100%, however, no discovery is possible. This observed impact of a 
large background uncertainty serves as an accuracy benchmark for background prediction 
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Sample 

20 % 

50% 

100 % 

0.1 

0.05 

0.04 

0.03 

0.1 

0.05 

0.04 

0.03 

0.1 

0.05 

0.04 

0.03 

SIMP (1 GeV) 

26 

> 

> 

> 

10 

32 

> 

> 

4.9 

16 

24 

34 

SIMP (10 GeV) 

26 

> 

> 

> 

10 

32 

> 

> 

4.9 

16 

24 

34 

SIMP (100 GeV) 

23 

> 

> 

> 

9.0 

29 

> 

> 

4.3 

14 

21 

30 

SIMP (200 GeV) 

16 

> 

> 

> 

6.4 

21 

31 

> 

2.9 

11 

16 

23 

SIMP (400 GeV) 

5.0 

20 

28 

36 

1.8 

8.4 

13 

18 

0.52 

4.0 

6.3 

9.4 

SIMP (700 GeV) 

0.33 

2.7 

4.6 

6.7 

< 

1.0 

2.0 

3.4 

<C 

0.11 

0.67 

1.5 

SIMP (1000 GeV) 

< 

0.13 

0.34 

0.64 

<C 

<C 

<C 

0.14 

<C 

<C 

<C 

<C 


TABLE IV. The Z-values of the signal samples with /'^qcd = 1 for a background uncertainty 
of 20%, 50% or 100% and CHEF cuts of 10%, 5%, 4% or 3%. The algorithm used to calculate 
these significances breaks down at very high or very low values, due to numerical imprecisions. 
This is represented by :§> and <C respectively. 


Sample 

20 % 

50% 

100 % 

0.1 

0.05 

0.04 

0.03 

0.1 

0.05 

0.04 

0.03 

0.1 

0.05 

0.04 

0.03 

SIMP (1 GeV) 

1.0 

6.1 

9.4 

13 

0.16 

2.5 

4.3 

6.7 

<c 

0.92 

1.9 

3.3 

SIMP (10 GeV) 

1.0 

6.1 

9.4 

13 

0.15 

2.5 

4.3 

6.6 

<c 

0.91 

1.9 

3.2 

SIMP (100 GeV) 

0.81 

5.1 

8.0 

11 

0.068 

2.1 

3.7 

5.8 

<c 

0.68 

1.6 

2.8 

SIMP (200 GeV) 

0.44 

3.3 

5.4 

7.9 

< 

1.3 

2.4 

4.0 

<c 

0.24 

0.89 

1.8 

SIMP (400 GeV) 

< 

0.73 

1.4 

2.2 

< 

0.077 

0.45 

1.0 

<c 

<C 

<C 

0.16 

SIMP (700 GeV) 

< 

< 

< 

0.057 

< 

<C 

< 

<C 

<c 

<c 

< 

< 

SIMP (1000 GeV) 

< 

< 

< 

<C 

< 

<c 

<C 

<C 

<c 

<c 

<C 

<C 


TABLE V. The Z-values of the signal samples with u^at/ctqcd = 0.1 for a background un¬ 
certainty of 20%, 50% or 100% and CHEF cuts of 10%, 5%, 4% or 3%. The algorithm used 
to calculate these significances breaks down at very high or very low values, due to numerical 
imprecisions. This is represented by ^ and <C respectively. 

methods in eventual data analysis. 


IV. ON COMPARISON WITH OTHER CONSTRAINTS 

Given the possibility of observing trackless jets at the LHC, it is tempting to try and 
set limits on the interactions of SIMPs and, possibly, to compare with other constraints. 
For this, we need to know 1) the high energy production cross section of SIMP pairs, 
as studied here, 2) the high energy inelastic cross of a SIMP with nucleus, , which 
is relevant for the response of the calorimeters (and some astrophysical constraints, like 
that from cosmic rays), and 3) the corresponding low energy elastic cross section, as 
reported in direct detection searches of dark matter, (and most of the astrophysical 
constraints). Given a simplihed SIMP model ([^, we may in principle consider calculating 
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these cross sections. This is provided the simplihed model holds over a very broad range 
of energies. Even if this is satished, we have to face the fact that the SIMP interactions 
involve large conplings, so that the validity of pertnrbative calcnlations is qnestionable. 
For those reasons, we consider the resnlts of the present section as tentative. Keeping this 
in mind, the main point we would like to explore is the following. Given that the SIMP is 
required to produce showers in the hadronic calorimeters, what is the low energy elastic 
cross section on nucleons? In absence of signal, this requirement sets a constraint on the 
SIMP cross section that may be compared with other searches. 

We begin with the problem of detection of the SIMP. The hadronic showers in the 
calorimeters will develop through inelastic scattering of the SIMPs with hadrons. As 
discussed in the previous section, our target value is a few mb, characteristic of standard 
hadronic cross sections. For instance the proton-proton total cross section at high energies 
(target frame) is about 40 mb, and roughly constant for proton momenta Plab ^ a 
few GeV [SB], while other hadrons have similar cross sections, in the range 10 — 40 mb 
(corresponding to a geometrical cross section with a radius of about 1 fermi), and similar 
energy dependence, exhibiting a plateau at high energies, plab GeV, which is the 
regime of deep inelastic scattering (DIS). A common feature is that, at these energies, the 
DIS cross section is always larger than the elastic one, which decreases with energy. At 
lower energies, the presence of resonances makes the situation much more complex and 
moreover dependent on the nature of the projectiles involved in the collisions, including 
neutrinos or, as a matter of fact, a SIMP. We will simply assume that the overall 
scale of inelastic scattering may be estimated from DIS; this should be conservative as 
the existence of resonances would increase the total cross section. In the case of a scalar 
mediator, we get from ([^ that the differential DIS cross section is given by 


d(j 

dxdy 


glgl uinE 
271 (Q2 + rniy 


Q^) + y'^xfiix, 


(scalar mediator) (4) 


where the functions /j (/j) are the parton distribution functions (PDFs) of the quarks 
(resp. antiquarks) inside the proton and x = /2m]sfi^, with u = E — E', and y = vjE 

the usual Bjorken variables, where E and E' are the energies of the incoming and scattered 
SIMPs in the lab frame. Similarly, for the vector mediator we have 


da^N 

dxdy 


g^gl m^E 
Att (Q^ -h m^)2 


^ xfi{x,Q‘^)+xfi{x,Q‘^){l-yf 


(vector mediator) (5) 


The behaviour of the total inelastic cross sections as a function of incoming energy E 
is similar for both mediators. They increase linearly with E for E < rri^/2mpf, and 
have a mild dependence on E for higher energies, being almost constant for the range of 
interest. This may be understood as follows. For high energies, large energy transfer is 
suppressed by the propagator, so that = 2m]\fXh' ~ m^, at which point the dominant 
contribution to the cross section comes from small x, x ^ m^/2mArF, with u E. A 
similar behaviour is predicted for ultra high energy neutrino DIS, in which case ~ 
p;0.4 [BB]- A major difference with neutrino DIS, is that for SIMPs the 

characteristic ~ may enter the non-perturbative regime if the mediator is very 
light, rUrj, < GeV: this is the main reason why we consider = 1 GeV. To the extent 
that the PDF may be trusted at small x and Q^j^from the expressions of the DIS cross 


® We have used the MSTW 2008 PDFs in this section |55] . 
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sections, we extract the couplings required for a given inelastic cross section. Taking for 
simplicity the couplings of the mediator to the light quarks and the DM to be equal, 

= Qg = g and requiring that = 10 mb for reference, we need ^ = 3 for = 1 
GeV for the vector mediator and g = 5.5 for = 1 GeV for the scalar mediator. We 
notice that, in the case of the scalar mediator, ( 7 ^/dvr ^ 1 , while they are still perturbative 
g‘^/4:7r < 1 for the vector mediator. The difference between the two cases stems from the 
different couplings to the PDFs. 

From this choice of parameters, we may now estimate the corresponding low energy 
elastic cross sections. From Q, the SIMP-nucleon elastic cross section through either the 
vector or the scalar mediator is given (in the Born approximation-see below) by 


C^xN - 


9 ^ 



( 6 ) 


where is the SIMP-nucleon reduced mass and is the effective coupling to the 
nucleons. For the scalar mediator, we have 


In = mN 




rrig 


+ — 


E 



(7) 


where f^g, dehned as m^f^g = {N\mgqq\N) , represents the contributions of the light 
quarks to the nucleon, and f'^Q = 1 — J2u,d,s frq [SD]- Using the values given in [HI], we 
get /tv ~ 11 (this is summing over all the quarks; taking only the hrst generation of 
quarks gives instead ~ 9.8). From the coupling required for inelastic scattering we 
get g^-N ~ 55, which is much larger than the acceptable value, g^^ ^ 6 , see section 
Hence the scalar mediator scenario is not viable, at least in the simplihed framework we 
consider. In the case of the vector mediator we have instead /jv = 3, which corresponds 
to g^N ~ 9. This is about 50 % larger than what we required, but the difference can be 
absorbed by adjusting the relative value of the DM and quarks couplings. 

This being said, we may question whether the Born approximation we used to calculate 
the low energy elastic cross section is valid for such large couplings. To check this, we have 
solved the Schrddinger equation for and yy scattering in the non-relativistic limit. 
This is akin to calculations of the so-called Sommerfeld effect, but here applied to repulsive 
interactions, see e.g. [2B]- Interestingly, the Sommerfeld effect is not only non-negligible 
but also helps in bringing the candidate in agreement with the astrophysical constraints on 
SIMP interaction with ordinary matter. The strongest limits have been reported in [3H], 
based on the impact of DM-baryon interactions on the matter power spectrum, which may 
be constrained using GMB anisotropies and Lyman-a measurements (see also 133 ). Our 
results for the elastic cross sections, and are shown in Figure]^ The hgure on 
the left shows the upper limit on from [58] together with the low energy cross section 
we have obtained solving the Schrddinger equation (the Born-approximation solution is 
shown for reference). The hgure on the right gives the corresponding self-interaction 
cross sectiorj^ together with the typical constraint from astrophysical observations, here 
from the Bullet cluster [HS]- It is important to appreciate that the cross sections are 
systematically smaller than their Born approximation value because our SIMP candidate 


® More rigorously, the so-called “transverse” cross section, as defined in [5^ . 
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Vector Mediator, Repulsive potential Mediator, Repulsive potential 



[GeV] [GeV] 


FIG. 3. Left) Elastic scattering for the vector mediator case with = 1 GeV with couplings 5 = 3. 
The perturbative (Born) approximation (blue line) is compared to non-perturbative result (green line). 
The short dashed line is the constraint from CMB+Lyman-ct of Ref.jSS]. Right) SIMP self-interaction for 
the same benchmark model of the left panel, comparing as well Born with the non-perturbative result. 
Colour code is identical to the left panel. The short-dashed line corresponds to an astrophysical reference 
limit on DM self-interactions cTxxI'^x ~ Icm^/g [SnilM]- Notice that the astrophysical constraint on 
interaction with baryons is more stringent than that on the self-interaction, as one would expect. 


has repulsive interactions]^ Notice that in both hgures we implicitly assume that the 
SIMP is the dominant form of DM. 

We now proceed to bring everything together. In Figure we show the relevant 
existing constraints in the a^N ~ '^x brief, the hgure shows the exclusion 

limits from direct detection searches, which are inoperant for large cross sections [21], 
those from various astrophysical observations that we just discussed and from rocket or 
balloon experiments (see Section 0 . and hnally those from DM searches at colliders. 
Most relevant for the present discussion is the line that corresponds to missing energy 
searches. The constraints are very strong for low mass DM candidates, complementary to 
those from direct searches. As the cross section of DM with nucleons increases however, at 
some point the signature is no longer missing energy, but trackless jets. This corresponds 
to the thick solid (red) line around a^N ~ 10“^® cm^ for > 1 GeV. If no excess above 
the QCD background is measured at the LHC, then the parameter space above this line 
would be excluded. The limit is conservative but robust, as a larger coupling implies a 
larger number of events (we took = 1 in Section HID), while the response of the detector 
would drop rapidly for a smaller coupling. We have tentatively extended the limit all the 
way up to m-x = 400 GeV, consistent with the analysis reported in Table IV and[V] while 
we are aware that our discussion of deep inelastic scattering is probably not completely 
reliable for such heavy SIMP candidates. As for the collider searches for WIMPs, the 
search of SIMPs is complementary with other constraints (essentially astrophysical in 
nature) for light candidates, below say 1 GeV. As is usual for such considerations, we 
emphasize again that the SIMP particles produced may not be the dominant form of DM, 
in which case the astrophysical constraints are essentially irrelevant. 


^ For self-interactions, the Born approximation breaks down if the particles are non-relativistic and 
5 ^/47r X > 1 [IH]. For scattering with nucleons, the dependence on rrix is replaced by the 

SIMP/nucleon reduced mass, rrix —t which explains why the non-perturbative result becomes 

independent of the SIMP mass for rux ^ 1 GeV). 
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FIG. 4. Summary plot showing all the most important applicable constraints. Our results are shown in 
the upper solid red line (“this work”), which corresponds to the green line of Figure (left). In black 
solid/dashed (lower lines), the monojet constraints are shown. The other constraints are: atmospheric 
XQC and RRS experiments (blue and cyan, respectively), underground experiments (brown dashed), and 
CMB-l-Lyman-a (black dashed). 


V. CONCLUSION 

In this work we have considered further the possibility that DM may be made (partially 
or totally) of particles with strong interactions with ordinary matter. These so-called 
SIMPs, for strongly interacting massive particles, are much less considered than their more 
popular siblings, the WIMPs, but they are regularly considered in the literature in order 
to address some astrophysical issues. While they are challenged by many observations, 
again mostly astrophysical, they are not completely excluded. Furthermore, little work 
has been done on possible constraints from colliders. Extending on previous works, in 
particular [11], we have studied in more details the possibility of observing trackless 
jets at the LHC, taking into account realistic simulations of the QCD background and 
the response of the detectors. Most notably, we show that the charged content of jets 
is a powerful discriminator to suppress dijet backgrounds at LHC, thus enhancing the 
sensitivity to a potential SIMP signal. Our analysis shows that SIMPs with mass up to 
rriy. ~ 400 GeV could lead to an observable signal, provided its interaction cross section 
with ordinary matter is about 10% of that of ordinary nucleons. Most of our work is 
dedicated to the forecast for the experimental search of SIMPs at the LHC. To do so, we 
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have adopted a simplified SIMP model. In this framework the SIMP interacts with quarks 
through a light mediator particle. This framework also allows, at least in principle, to 
compare the sensitivity reach of the LHC search with other constraints, most of which 
are at much lower energies. We show this in Figure]^ assuming that our SIMP simplified 
candidate constitutes the dominant form of dark matter. Where relevant, we pointed out 
the limitations due to the difficulty of doing reliable calculations with strongly interacting 
particles, be them within a simplified framework. The signature and the proposed analysis 
are however essentially independent of these potential complications, i.e. they stand 
by themselves. While much improvement could also be envisioned for the experimental 
analysis, most notably a more realistic analysis of the detector response, our feasibility 
study shows that SIMP candidates over a wide mass range may be efficiently searched at 
the LHC. 
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